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Abstract 

Many chaotic dynamical systems of physical interest present a strong form of non-hyperbolicity 
called unstable dimension variability (UDV), for which the chaotic invariant set contains periodic 
orbits possessing different number of unstable eigendirections. The onset of UDV is usually related 
to the loss of transversal stability of an unstable fixed point embedded in the chaotic set. In this 
paper, we present a new mechanism for the onset of UDV, whereby the period of the unstable 
orbits losing transversal stability tends to infinity as we approach the onset of UDV. This mech- 
anism is unveiled by means of a periodic orbit analysis of the invariant chaotic attractor for two 
model dynamical systems with phase spaces of low dimensionality, and seems to depend heavily 
on the chaotic dynamics in the invariant set. We also described, for these systems, the blowout 
bifurcation (for which the chaotic set as a whole loses transversal stability) and its relation with 
the situation where the effects of UDV are the most intense. For the latter point, we found that 
chaotic trajectories off but very close to the invariant set exhibit the same scaling characteristic of 
the so-called on-off intermittency. 
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I. INTRODUCTION 



A chaotic trajectory of a dynamical system is structured upon an infinite set 
of unstable periodic orbits. A systematic investigation of the latter would pro- 
vide an understanding of the ergodic properties of the chaotic attractor itself, 
like its Lyapunov exponents and multifractal spectrum. The relation between 
chaotic trajectories and unstable periodic orbits embedded therein has been ex- 
plored in some detail for hyperbolic dynamical systems and non-hyperbolic ones 
when there are homoclinic tangencies. There is, however, a stronger form of 
non-hyperbolicity called unstable dimension variability (UDV), which is charac- 
terized by the coexistence, in the chaotic invariant set, of periodic orbits with a 
different number of unstable directions. We harnessed the connection between 
unstable orbits and chaotic trajectories to investigate the parametric evolution 
of a chaotic set presenting unstable dimension variability. The onset of UDV 
is usually related to the loss of transversal stability of an unstable low-period 
saddle embedded in the chaotic set, through a local bifurcation. We have found, 
however, a new and qualitatively different scenario underlying the onset of UDV: 
as we approach this point, the period of the orbits embedded in the chaotic set 
and losing transversal stability tends to infinity. This fact seems to depend in 
a sensitive way on the chaotic dynamics in the inveriant set. Another point of 
interest is the blowout bifurcation point, where the chaotic set as a whole loses 
transversal stability, the effects of UDV being the most severe in terms of the 
shadowing properties of computer-generated trajectories. In general, once we 
have UDV the shadowing distances between fiducial chaotic orbits and noisy 
ones diverge very fast (i.e. the shadowing times may become small). We made 
a periodic orbit analysis of the chaotic attractors for two low-dimensional dis- 
crete dynamical systems as paradigmatic models, where this new mechanism 
for the onset of UDV, as well as the blowout bifurcation point, are described. 
We expect that UDV be rather common in high-dimensional dynamical systems 
of physical and biological interest, like lattices of coupled oscillators or maps, 
hence it is important to characterize observable manifestations of the existence 
of UDV in a given system. Three of them are treated in this paper, in order to 
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provide numerical fingerprints of the scenarios for UDV we have described: (i) 
the statistics of finite-time Lyapunov exponents (in the transversal direction); 
(ii) the intermittency exhibited by chaotic orbits off but very close the chaotic 
set; and (iii) the distribution of log-shadowing distances to the chaotic set. In 
the first place, at the onset of UDV, the fraction of positive finite-time transver- 
sal exponents becomes nonzero, while it takes on the value 1/2 for the blowout 
bifurcation point. Secondly, the intermittent switching between a quiescent and 
a bursting state is influenced by the transversely unstable orbits embedded in 
the chaotic set, and it has the universal features of the on-off intermittency as 
we approach the blowout bifurcation point. Finally, we characterized statisti- 
cally UDV in terms of the exponential distribution of log-(shadowing) distances 
to the invariant set, and the corresponding power-law scaling exhibited by the 
corresponding shadowing times. 

Unstable periodic orbits (UPOs) are the soul of dynamics, as has been already noted by 
Poincare in his work on the restricted three-body problem [1]. More specifically, UPOs form 
the skeleton upon which the chaotic dynamics is built. In fact, there are infinitely many 
UPOs embedded in a chaotic attractor. From the ergodic point of view, this description 
amounts to relate the natural measure generated by a typical chaotic trajectory to the 
infinite number of atypical measures generated by each of the UPOs embedded in the chaotic 
attractor [2], This connection has been rigorously demonstrated for hyperbolic dynamical 
systems, for which the splitting into stable and unstable invariant subspaces is consistent 
under the evolution of the dynamics [3]. But it has also been numerically verified for two- 
dimensional non-hyperbolic systems possessing homoclinic tangencies, which violate the 
invariant space splitting at the tangency points [4], 

There is another class of non-hyperbolic dynamical systems, whose invariant sets also fail 
to exhibit a consistent invariant subspace splitting. As the trajectory evolves on the chaotic 
attractor, the dimension of the unstable subspace (and hence the stable) varies densely [5]. 
In this situation, the breakdown of hyperbolicity is even more severe than when homoclinic 
tangencies are present, since in the latter case we may found shadowing trajectories for a 
time often large enough to give meaning to computer-generated single chaotic trajectories 
[6]. When the unstable dimension varies, however, shadowing trajectories may exist for 
a time rather short so that in practice one cannot be sure of the validity of the chaotic 
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trajectories numerically generated [7]. For these systems, whose invariant chaotic set is 
said to exhibit unstable dimension variability (UDV), the general relation between natural 
measure and UPOs is still an open problem. 

Two key problems arise in the study of UDV in dynamical systems: the first is its 
onset, normally related to a bifurcation where some low-period UPO (and its infinitely 
numerous pre-images) change its unstable dimension [8]. The second problem addresses 
the point where UDV is the most intense, usually referred to as a blowout bifurcation, 
since it marks the point at which a chaotic attractor lying in an invariant subspace loses 
transversal stability as a whole [9-12]. However, not only systems possessing such invariant 
subspace can exhibit UDV. One counter-example is the kicked double rotor map [13], which 
was found to exhibit UDV for a limited parameter range, in spite of not having such an 
invariant manifold embedded in its four-dimensional phase space [14]. One of the numerical 
fingerprints of UDV is the fluctuations about zero observed for the finite-time Lyapunov 
exponent (FTLE) closest to zero [15]. 

In this paper we study the loss of transversal stability of the chaotic invariant set and its 
relation with the UPOs embedded in it, interpreting our results in terms of the hyperbolicity 
breakdown caused by UDV. Two archetypical dynamical systems are used to investigate the 
dynamics: a two-dimensional and a three-dimensional map, possessing invariant subspaces 
which are chaotic. The blowout bifurcation is interpreted in terms of the parameter value 
for which the effect of UDV is the most intense [16]. On the other hand, the onset of UDV 
is described by a mechanism hitherto not considered before: as we approach this point the 
period of the orbits losing transversal stability tends to infinity. In previous studies [17, 18] 
the onset of UDV has been described through the loss of transversal stability of a low-period 
orbit embedded in the chaotic set. We found that the particular scenario for the onset of 
UDV depends on the chaotic dynamics in the invariant set. 

The contribution of the UPOs embedded in the chaotic attractor is numerically analyzed 
by comparing the fraction of positive transversal finite-time Lyapunov exponents (FTLE) 
for typical chaotic orbits with the fraction of transversely unstable periodic orbits embedded 
in the chaotic set [16]. For both system considered in this paper, we find that the onset of 
UDV and the blowout bifurcation correspond to those points for which this fraction becomes 
nonzero and it is equal to 1/2, respectively [19]. 

We also consider the dynamics of chaotic trajectories off but very near the invariant 
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chaotic set. It turns out that they are strongly influenced by the transversal stability of 
the UPOs embedded in the chaotic set. One manifestation of this influence is a UDV- 
induced intermittency: the trajectories remain close to the invariant subspace until they 
approach a transversely unstable orbit and burst away for a while before returning to the 
vicinity of the invariant subspace [20]. We propose an interpretation of this phenomenon 
from the point of view of shadowing of chaotic trajectories in the invariant set due to UDV: 
the shadowing times are the lengths of the laminar intervals for which trajectories remain 
close to the invariant subspace [21]. The shadowing log-distances and times are found to 
have exponencial and power-law probability distributions, respectively, in accordance with 
a biased random- walk model with reflecting barrier [7]. 

The remaining of this paper is organized as follows: Section II introduces the systems we 
use to perform the numerical experiments, whereas Section III deals with the computation 
of the corresponding finite-time Lyapunov exponents and how they can be used to describe 
UDV. Section IV presents a periodic orbit description of the chaotic attractors focusing on 
the measures generated by UPOs with different unstable dimensions. Section V discusses 
the UDV-induced intermittency observed for chaotic trajectories off the invariant set and the 
statistical properties of log-shadowing distances and shadowing times for such trajectories. 
Our conclusions are left for the last Section. 



II. DYNAMICAL SYSTEMS AND THEIR PROPERTIES 

There are extensive numerical evidence that non-hyperbolic systems with UDV appear 
typically in prototype models of spatially extended dynamical systems, like coupled maps 
or oscillator lattices [22]. However, due to the large dimensionality of their phase spaces, 
a periodic-orbit treatment may become a quite time-consuming task. Hence, we choose to 
analyze simpler dynamical systems but which retain most of the essential dynamical features 
expected from those complex lattices. 

A typical TV-dimensional system with such characteristics is [4, 11, 23] 

x n+ i = f(x n ), (1) 
z n+ i = F(x„ ) a)G(z n ) ) (2) 

where x G R a 'ii (N\\ > 1), and z G R N± (N ± > 1) such that N = N\\ + 7V_l, and a > 
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is a bifurcation parameter. We assume the existence of a symmetry in the system so as to 
warrant the condition G(— z) = — G(z) and, as a consequence, that G(0) = 0, which we 
interpret as defining a iVy-dimensional invariant subspace S given by z = 0. In terms of a 
coupled lattice map, for example, S can stand for the synchronization manifold wherein the 
time evolution is the same for each map. Our goal of describing the parametric evolution of 
UDV demands that the dynamics in S be governed by a map f (x n ) with a chaotic attractor. 
The jV_L-dimensional vector z n represents the directions transversal to the invariant manifold, 
which turn to be important to characterize the stability of the completely synchronized state. 

In this paper we analyze two different possibilities for the chaotic dynamics in the invari- 
ant subspace [4, 23]: 

• Model B, where JVy = 1, for which 

x n+1 = 2x n (mod 1), (3) 

or a Bernoulli shift with a dense chaotic orbit in [0, 1) with Lyapunov exponent X x = 
In 2. In order to generate the numerical trajectories the map was slightly modified to 
x n+ i = 2x n + (35 x mod 1, where (3 is or 1 with equal probabilities and 5 X = 2 -48 . Due 
to the transitive property of the strong chaos in [0, 1), the trajectories so generated, 
although noisy, are shadowed by true chaotic trajectories of the Bernoulli shift map 
for an arbitrarily long interval [24]. 

• Model C, where N\\ = 2, for which 

= 2x n + y n (mod 1), (4) 

= Xn + Vn (mod 1), (5) 

is the Arnol'd cat map with a hyperbolic chaotic set in the unit square [0, 1) x [0, 1) 
[25]. 

We will consider a one-dimensional transversal dynamics (iVj_ = 1) given, for both models, 
by the following realization of Eq. (2) 

= ax n g(z n ), (6) 

where we impose that g(0) = and g'(0) = const. Since x n has a chaotic behavior, the 
function F(x n ] a) = ax n stands for a driving from the x-dynamics in the invariant subspace 
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to the symmetric z-subsystem. We have used the function g(z) = (1/2tt) sin(27rz), but other 
maps can be used as well, like the logistic map g(z) = z(l — z). 

The Jacobian matrices for models B and C are given, respectively, by 

(2 \ 

DB(x„) = , (7) 

\ (a/27r) sin(27T2: n ) ax n cos(27rz n ) / 

and 

^ 2 10^ 

DC(x n )= 110. (8) 

^ (a/27r) sin(27rz n ) ax n cos(27rz n ) j 
The Bernoulli shift map (describing the dynamics in the invariant manifold of map B) 
has N p = TP — 1 unstable periodic points of period p given by 

ZpO") = 2?3i' (j = 0,l > ...2'-2). (9) 

Each period-p orbit is represented by x p (k), where k is such that x p (k) = mf{B l (x p (k))} 3 ?~Q. 
Using this notation the fixed point is Xi(0) = 0. When p = 5, for example, there are 
out of 2 5 - 1 = 31 period-5 points, for example x 5 (l) = {1/31,2/31,4/31,8/31,16/31}, 
£5(11) = {13/31, 26/31, 21/31, 11/31, 22/31}, and so on. 

Similarly, for the cat map, which describes the dynamics in the invariant manifold of C, 
the counting of the periodic orbits involves F p , the pth term of the Fibonacci sequence, given 
by Fq = F\ = 1 and F p = F p _i + F^_ 2 (for p > 2). The number of period-p (unstable) points 
is [26] 

N p = (F p + F p _ 2 f - 2[1 + (-1)% (10) 
their values being solutions of 

Zp(J) = A p x pU) + BpVpti) (modi), (11) 
y P (j) = C p x p (j) + D p y p (j) (mod 1), (12) 

where A p = F 2p , B p — C p = F 2p -\, D p = F 2p -2- These solutions are non-negative rational 
numbers, for which x p (j) = a/c and y p (j) = b/c, for some a, b = 0, 1, 2, ... c — 1 and we have 
defined [26] 

(F p + F p _ 2 , ifp>lisodd, 
(13) 

y5F p _i, if p > 1 is even. 
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The reason for choosing dynamical systems presenting an invariant subspace in which 
there is a hyperbolic chaotic set will be clarified in the next sections, where we aim to relate 
the natural measure of this set with the expanding eigenvalues of the unstable periodic orbits 
embedded in the chaotic set. 

III. FINITE TIME LYAPUNOV EXPONENTS 

The dynamical systems introduced in the last section have some common points of interest 
that are relevant for our work. Both are skew-symmetric and have invariant subspaces in 
which the dynamics is strongly chaotic; they are simple enough to allow for a counting 
of unstable periodic orbits, and they both exhibit the transition to unstable dimension 
variability through a bifurcation of some periodic orbit. Such a bifurcation changes the 
unstable dimension of the periodic orbit and its pre-images. As a result of this bifurcation, 
in the neighborhood of the invariant manifold there appear tongue-like structures anchored 
at the invariant subspace. The open dense set of tongues are associated with the particular 
periodic orbit that underwent the bifurcation. The closed set is the stable manifold of the 
remaining periodic orbits and of the chaotic orbits. 

A way to quantify the relative abundance of periodic orbits with a different number of 
unstable directions is to compute the finite, or time-n Lyapunov exponents of typical chaotic 
orbits. They are defined, for a D-dimensional map F(x) [27], 

Ai(xo,rO = -ln||DF"(x ).e i || (i = 1,2,...L>), (14) 

n ....... 

where is the zth singular vector corresponding to the singular value £™ of the Jacobian 
matrix DF n of the nth iterated map. The usual Lyapunov exponents are obtained as the 
infinite-time limit of the above expression. While, in general, the time-ra exponent depends 
on the initial condition, its infinite-time limit, does not, except for a set of points with 
Lebesgue measure zero: 

Xi = lim Ai(x 0j n). (15) 

The eigenvalues of the product of the Jacobian matrices, evaluated at a given point, 
DB™(x ), are given by 

n-l 

£ = 2", e = a"II^ C0S ( 27r ^)' ( 16 ) 

i=0 

8 



and, for DC™(xq), are 




(17) 



n-l 




i=0 



The characterization of non-hyperbolicity via unstable dimension variability is based on 
the finite-time exponent closest to zero. For our low-dimensional maps, this turns to be 
the exponent computed along the transversal direction to the invariant manifold, which 
we henceforth will refer to as \±(n). The time-n transversal exponent for both models, 
according to Eqs. (16) and (17), is given by 



such that the transversal exponent is the infinite-time limit of the above expression: 



where we have used the Poincare ergo die theorem and p(x) = 1 is the invariant density 
of iterations in the unit interval for model B or its two-dimensional equivalent for the unit 
square in model C. It follows from Eq. (19) that the blowout bifurcation, i.e., the point 
where the invariant manifold becomes transversely unstable as a whole, occurs at a* = e for 
both models. 

A fingerprint of unstable dimension variability (UDV) is the existence of positive and 
negative fluctuations for X±(n). This means that the transversal exponents along a chaotic 
trajectory typically present finite-time segments for which the phase-space volumes are lo- 
cally either expanding or shrinking on average. We take a long chaotic trajectory and 
partition it into a large number of segments with length n. After computing the time-n ex- 
ponents for each trajectory segment, we are able to obtain a frequency distribution, related 
to a probability density P(X_\_(n)). For n large enough this distribution has a Gaussian-like 
shape, whose width decreases as ~ n~ 1 / 2 such that, in the limit n — > oo it reduces to a 
delta- function centered at the infinite-limit A_l [19]. It should be mentioned, however, that 
UDV is not necessarily related with the non-Gaussian nature of the probability distributions 
but with the fact that they have tails of positive finite-time exponents. 




n-l 



(18) 



i=0 




(19) 
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FIG. 1: Probability distribution of time-n transversal Lyapunov exponents for maps (a) B, n — 28; 
and (b) C, n = 16, and different values of the bifurcation parameter a. Solid lines: distribu- 
tions evaluated using typical chaotic trajectories with n x N n points. Dashed lines: distributions 
computed using all unstable period-n orbits. 
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TABLE I: Model B 

& ^max X± 

1.00 - 0.86 -1.00 

2.72 0.14 0.00 

6.00 0.79 0.92 



TABLE II: Model C 

OL Xmax A_L 

1.00 - 0.99 -0.96 

2.72 0.00 0.06 

6.00 0.76 0.84 



Numerical approximations to this probability density for time-ra transversal exponents are 
shown as solid lines in Figures 1(a) and (b) for models B (n = 28) and C (n = 16), respec- 
tively. The dashed line represents the analogous distribution for all UPOs of period-p (p = n) 
embedded in the chaotic attractor, which nearly coincides with the distributions obtained 
from the typical chaotic trajectories. We found both distributions to have a Gaussian-like 
shape, drifting without noticeable distortion towards positive exponents as the control pa- 
rameter a increases [28]. Actually the distributions have uneven tails, the positive tail having 
a cutoff at the point X±(n) = In a, whereas the negative tail looks like a Gaussian one. If 
the distributions were Gaussians their maxima A max would coincide with the infinite-time 
exponent A_l given by Eq. (19). This is not the case for the distributions shown in Fig. 1, 
though, as can be seen in Tables I and II where we present the corresponding values of X max 
and A_l (obtained from Eq. (19), which are clearly different. 

There are two critical points to be explored in the evolution of the distribution P(X±(n)) 
for increasing a: (i) if a < a c & 1.0 there are no positive fluctuations of the finite-time 
exponent, indicating absence of UDV. The critical value a c signals the onset of UDV, since 
positive fluctuations start appearing at this point due to the existence of transversely un- 
stable UPOs embedded in the attractor; (ii) for a c < a < a* = e, there are more negative 
than positive fluctuations, indicating that the distribution peak lies in the negative semi- 
axis. Accordingly, before the blowout bifurcation the infinite-time exponent A_l is likewise 
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negative, but, at the blowout bifurcation a* = e, the distribution peak is slightly displaced 
from zero, since the distributions themselves have unequal tails due to non-gaussianity. (hi) 
When a > a* the positive fluctuations dominate. 

Another way to understand these observations is to compute the positive fraction of the 
transversal finite-time exponent, or 



Mn) = ^— V V - — 20 

The solid lines depicted in Figures 2(a) and (b) present the computed fraction of time-ra 
transversal exponents for maps B (n = 28) and C (n = 16), respectively, as a function of the 
bifurcation parameter a. Since this fraction exhibits a monotonically increasing trend with 
a, we can use the condition <fi > for characterizing the onset of UDV in those systems, 
giving an estimated value of a c « 1.0 for this critical point. However, the insets in Fig. 2 
reveal small fluctuations near this critical point (due to the finite bin width) making difficult 
a precise evaluation of a c based solely on this criterion. This is not surprising, however, since 
the characterization of UDV through finite-time Lyapunov exponents is difficult, inasmuch, 
we must deal with exponents fluctuating about zero, and at this parameter value the density 
of the attractor is small so that very long trajectories are needed to obtain reliable results. In 
order to overcome such problems, a numerical method based on the scaling of the probability 
distribution of points near the fixed points has been recently proposed [29]. 

At the blowout bifurcation, there follows from Fig. 2 that nearly half of the finite-time 
transversal exponents are positive: this fraction actually takes on values around 0.57 and 0.54 
for models B and C, respectively. The difference with the theoretical value 0.5 is due to the 
use of a finite (and small) time interval n in the computation of the probability distribution 
X±(n). We remark that we have used such small values of n since we set n equal to the period 
p of the UPO's used as examples of atypical orbits. We are particularly interested to explore 
the situation that lies between the onset of UDV, for which the first UPO loses transversal 
stability, and the blowout bifurcation, where the manifold as a whole becomes transversely 
unstable. However, this does not necessarily mean that for the blowout bifurcation all UPOs 
have lost transversal stability. The answer lies in the relative contribution of the two sets of 
UPOs with different unstable dimensions. 
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FIG. 2: Fraction of positive time-n transversal Lyapunov exponents for maps (a) B (n = 28) and 
(b) C (n = 16) as a function of the bifurcation parameter a. Solid lines: typical chaotic orbits; 
dashed lines: fraction of transversely unstable period-p orbits. The insets in both figures show 
blowouts of the regions near the onset of UDV for both maps. 
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IV. TRANSVERSAL STABILITY OF UNSTABLE PERIODIC ORBITS 



After the onset of UDV (a > a c ) there are two kinds of coexisting UPOs: saddles and 
repellers. The saddles (repellers) are transversely stable (unstable) and thus have unstable 
dimension one (two). Just after the onset of UDV, the contribution of repellers to the 
measure of the chaotic set of the invariant subspace is very weak, since the UPO which lost 
transversal stability has Lebesgue measure zero. As one approaches the blowout bifurcation 
at a* the contribution of repellers becomes increasingly significant until it matches the 
contribution of the saddles. After the blowout bifurcation point, the contribution of repellers 
to the measure become dominant. 

The relative contributions from saddles and repellers in the context of UDV can be 
made more precise in the context of periodic-orbit theory [11, 12]. In order to compute the 
infinite-time transverse Lyapunov exponent A_l, we use typical trajectories on the chaotic 
invariant set with respect to its natural measure ^(^) (^ is the unit interval [0, 1] and 
the unit square for models B and C, respectively). Since there are an infinite number of 
UPOs embedded in fi, they support the natural measure in the sense that, when computing 
A_l, these periodic orbits contribute with different weights. These weights are given by the 
natural measure of a typical trajectory which visits the neighborhoods of the periodic orbits, 
and are related to the magnitude of the unstable eigenvalues of those unstable orbits [3]. 

The influence of the UPO's on the loss of transversal stability of the chaotic invariant 
set can be appreciated in Fig. 3(a), where we plotted the critical value of the bifurcation 
parameter af irst (p) for which the first UPO of period p loses transversal stability, as a 
function of the period p of the UPO. The value of af irst (p) decreases monotonically for 
increasing period and seems to asymptote to 1 for large p, this evolution being essentially 
the same for models B and C; a result which is in close agreement with the estimate we 
made from Fig. 2. These observations can be made rigorous by a straightforward stability 
calculation we present in detail in an Appendix. We conclude that the first orbits to lose 
transversal stability are those with very large periods. In this sense, the atypical trajectories 
represented by UPOs do not dominate the effects produced by typical trajectories in the 
chaotic invariant set. 

Figure 3(b) depicts the values of a for which the last UPO of period-p loses transversal 
stability for models B (circles) and C (crosses). The solid lines are least squares fits of the 
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form exp(^p), where tp = 0.35 and 0.26 for models B and C, respectively. If we consider a 
fixed value of a, say a = 2, it would seem that most UPOs (mainly with large period) have 
already lost transversal stability. However, Fig. 3(b) shows that, at a = 2, there are still 
transversely stable UPOs of all periods. The conclusion we get is that large-period UPOs 
lose transversal stability earlier (in terms of increasing a-values) than low-period ones, but 
it takes longer to exhaust all possibilities. 

This can be viewed as a consequence of the number of period-p orbits to increase with 
their period p as e hTP , where hx > is the topological entropy of the chaotic set, such 
that there are far more UPOs with large periods than low-period ones [30]. A standard 
calculation shows that = In 2 ^ 0.69 for model B, and 



for model C. Hence the number of period-p orbits increases more rapidly for model C than 
for model B, such that, for a fixed p, it is more difficult to count these orbits for model C. For 
this reason, we have limited ourselves to UPOs of periods up to 28 and 16 when considering 
models B and C, respectively. Moreover, we also expect predictions from periodic orbit 
theory to agree more with numerical results for model B. 

The fact that the onset of UDV is marked by the loss of transversal stability of a UPO 
with infinite period seems to depend in a sensitive way on the chaotic dynamics in the 
invariant set, keeping the same transversal dynamics. For example, when we used in model 
B, a tent map x n ^i = 1 — 2\x n — (1/2) | (or its topological conjugate, = 4x n (l — x n )) 

instead of the Bernoulli shift map, the first orbit to lose transversal stability is a fixed point. 

The natural measure of a typical trajectory in the neighborhood of a periodic orbit is 
related to the probability of being in its close vicinity, which is smaller the more unstable 
the periodic orbit is [2]. Hence UPOs with large unstable eigenvalues have a comparatively 
smaller contribution to the natural measure. Summing over all unstable period-p orbits 
embedded in the invariant set Q gives then its natural measure [3] 



where x p (j) is the jth fixed point of F p (x) belonging to the set Q, i.e., x p (j) is on a period-r 
orbit, where r equals either p or a prime factor of p. 




(21) 




(22) 
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p 



FIG. 3: Critical values a p of the bifurcation parameter for which the (a) first and (b) last period-p 
orbits lose transversal stability. Circles and crosses refer to models B and C, respectively. The lines 
in (b) are least squares fits with slopes 0.35 and 0.26 for models B and C, respectively. 
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According to Ref. [3], L u is the absolute value of the product of eigenvalues related to 
the unstable directions of DF P , computed for the orbit points x p (j): 

D 



i=l 



(23) 



such that > 1 and Z? is the dimensionality of the map F. For the models treated in 
this paper, it turns out that L u does not depend on x p (j) since the x-direction is always 
unstable and £ x takes on the same values for all orbits x p (j), such that L u = 2 P for model 
B, and L u = (1/2) (F 2p + F 2p _ 1 + V5F 2p _ 2 ) for model C. 

Equation (22) was originally derived for hyperbolic chaotic sets, for which there are 
Markov partitions formed by intersections of stable and unstable manifolds of UPO embed- 
ded in the chaotic invariant set Q [3]. In subsequent papers [4, 11] non-hyperbolic systems 
with homo clinic tangencies between stable and unstable manifolds of unstable periodic or- 
bits were considered, and Eq. (22) was used by assuming that L u would be the absolute 
value of the greater eigenvalue of the orbit x p (j), or, for model C, 

L w (x p (j)) = su P {ia Kjugi}. (24) 

In many instances, Eqs. (23) and (24) would give the same result for the natural measure 
of typical chaotic orbits, ^(fi), provided one computes only the eigenvalues related to unstable 
directions within the invariant set Q. The results of Refs. [4, 11] did not differ from the 
original definition of L U) Eq. (23), because there was only one unstable direction for orbits 
in the chaotic set Q. On the other hand, higher dimensional systems such as the ones we 
deal with in this paper, there can be unstable directions not in the chaotic invariant set 

namely those related to transversal directions. Since these chaotic sets are taken to be 
hyperbolic by construction, Eq. (22) can be applied to them, but there may be differences 
between the results furnished by Eqs. (23) and (24). However, since an unstable direction 
(the x-direction) along the invariant subspace has usually the greater eigenvalue, the possible 
discrepancies between (23) and (24) are conceptually significant though small in size as we 
shall show. 

The (atypical) natural measure associated with the jth period-p orbit is the normalized 
ratio [16] 

= jV 1/L " (X ^°" )) , (25) 
17 



where N p is the number of period-p orbits. N* and are the numbers of transversely 
stable and unstable period-g orbits, respectively, such that N* + = N p . The weights of 
the transversely stable and unstable period-g orbits are given, respectively, by 

A; = ^^0')A±(x p (i),p) (for A ± (x P 0%p)<0), (26) 
A; = J]^0')A^(x p (j),p) (for A ± (x P 0') ) p)>0), (27) 

where \j_(x p (j),p) is the time-p transversal Lyapunov exponent for the jth period-p orbit. 
If \±{x p (j),p) is positive (negative) the periodic orbit is transversely unstable (stable). 

Within the realm of the periodic-orbit theory, when A_l = (blowout bifurcation), it fol- 
lows that the contributions of the transversely stable and unstable period-g orbits are exactly 
matched. In this sense, we can say that UDV is the most intense at blowout bifurcation, and 
the validity of computer-generated chaotic trajectories is rather short since there would be 
no shadowing trajectories of the chaotic invariant set for a significant lifespan [31]. In terms 
of the finite-time transverse exponents, one could expect that half of the time-ra exponents 
would be positive, which would imply that, at the blowout bifurcation, <p(n) = 1/2. While 
this connection has been numerically verified for some dynamical systems [19], an inspection 
of Fig. 2 reveals discrepancies with this assertion, for the blowout bifurcation occurs slightly 
after <fi equals 0.5. 

Another way of analyzing the relative contribution of period-g saddles and repellers is to 
use the so-called contrast measure, defined as 

C p = |/£-a£|, ( 28 ) 

where 

^ = X>p«, /x; = X> P (0 = l-A^, (29) 

i=l i=l 

and fi p (i) is given by Eq. (25). The variation of the contrast measure with a is depicted in 
Fig. 4 for UPOs with period-p. Before the onset of UDV (a < a c ) the only contribution is 
that from the saddles, hence C p = 1. The vanishing of C p marks the point where the relative 
contributions of period-p saddles and repellers exactly match each other. For both models, 
this point does not coincide with the blowout bifurcation at a* = e, although for model C the 
difference is smaller than for model B, in agreement with the results obtained for the fraction 
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FIG. 4: Contrast measure versus bifurcation parameter for unstable orbits with period: (a) p = 28 
(model B); (b) p = 16 (model C). 

of positive finite time exponents in Fig. 2. The explanation for this small discrepancy lies 
in the use we made of finite and small periods for the UPOs to which atypical measures are 
compared. We expect periodic-orbit predictions to match better with numerical results if p 
is enlarged. 

Although analytical methods to count unstable periodic orbits exist only for low- 
dimensional maps, such as those studied in this paper, there are powerful and efficient 
numerical methods devised to do this task even for more complex systems [32-34] The prob- 
lem with those methods is that the computational time required grows both with the period 
of the orbit as well as with the dimension of the dynamical system. For a D-dimensional 
map, the method proposed in Ref. [34] requires that, in order to find orbits of period p+ 1, 
out of 2 du recursive equations (involving D x D matrices) be solved for each period-p orbit 
with unstable dimension d u . But the number of periodic points grows exponentially fast 
according to the topological entropy. As a representative example of high-dimensional dy- 
namical system, for the kicked double rotor map [13] there are more than 300, 000 period-7 
points with unstable dimensions equal to 1 or 2, such that one must deal with approximately 
a million of matrix equations to find the period-8 orbits. On the other hand, since in the 
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FIG. 5: Average number of iterations for a typical trajectory of model B (circles) and C (crosses) 
to be trapped in a vicinity of width 10 -12 of the chaotic set versus the difference between the 
bifurcation parameter and its critical value a* = e. The solid lines are least squares fits, cf. Eq. 
(30), with slopes -0.852 and -0.447 for model B, and -0.885 and -0.177 for model C. 

scenario immediately following the onset of UDV we need orbits of arbitrarily large period, 
the use of such numerical methods might not be feasible to find the value of a c . That was 
essentially the reason behind our using of simple, low-dimensional maps for which a direct 
analytical treatment is feasible. In order to establish the existence of such scenario in more 
complex systems, like higher-dimensional flows or coupled maps, one would need to resort 
to those time-consuming methods. 

V. UDV-INDUCED INTERMITTENCY 

The existence of UDV in the invariant chaotic set Q for the models considered in this 
work, has an observable effect on trajectories which do not belong exactly to the chaotic 
set but are near enough to feel the coexistence of periodic orbits transversely stable and 
unstable. These trajectories have typically an intermittent behavior, alternating quiescent 
and bursting phases. We can call the process UDV-induced intermittency, since here chaotic 




20 



bursting is accompanied by the lack of hyperbolicity [20]. 

A trajectory which is initialized off but very close to the invariant chaotic set Q will 
be pushed either away or towards Q if it is in the vicinity of a repeller or a saddle UPO, 
respectively. The time it takes for such a trajectory to leave a small neighborhood of H, say, 
of size 10~ 12 , diverges rapidly as we are close to the blowout bifurcation, as shown by Fig. 
5, where we observe two scaling regions for which 

<r >~ \a-a*\~™ (30) 

where the slopes are different for the region near and far the critical point, with a shoulder 
at \a — a* | = 0(1). As a consequence, if we are interested to investigate the behavior of 
long trajectories in the vicinity of blowout bifurcation, it would take a prohibitively long 
time unless we add an extrinsic noise ^R n in the transversal dynamics of both models B and 
C, where 7 is a noise level of 10~ 16 and R n is a pseudo-random variable in the interval [0, 1] 
[7, 21]. 

This procedure does not alter the transversal dynamics in the scenario of UDV, since 
we can regard the transversal finite-time Lyapunov exponents as random kicks that push 
the trajectory away or towards Q as the exponents are positive or negative, respectively. 
Hence we are just adding an extra noisy term to an already noisy system (the transversal 
dynamics near a*, where the positive and negative finite-time Lyapunov exponents match 
each other) that does not alter the difference between up and down kicks due to UDV. 
However, this procedure should not be used too far from the blowout bifurcation, otherwise 
the distribution of kicks due to the finite-time exponents is biased towards either positive or 
negative value, and adding a white noise with uniform distribution would lead to artificial 
results. 

Representative examples of intermittent behavior induced by UDV are depicted in Figs. 
6(a) and (b) for models B and C, respectively, and a value of a very close to the blowout 
bifurcation: a = a* + <5, with 5 ^ 10~ 3 . In both cases we can define a laminar phase as the 
inter-burst interval with duration r. As in most intermittent scenarios, the duration of these 
laminar phases differ greatly but they are statistically distributed according to a power-law 
scaling [cf. Fig. 7, where the scaling was obtained through a least-squares fit] P(r) ~ 
where u was estimated to take on the values 1.492 and 1.480 for maps B and C, respectively. 
Both values compare well with the theoretical value 3/2 predicted for on-off intermittency 
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FIG. 6: Evolution of the transversal coordinate of a trajectory off but very close to the chaotic set 
of models (a) B and (b) C, with a = 2.72 > a*. 

[35]. 

The connection between UDV and on-off intermittency has been pointed out in Refs. [36], 
where it is claimed to be a universal feature. In fact, the observed intermittent transition 
to complete synchronization of coupled maps satisfies the same power-law scaling with 3/2 
exponent [18, 31]. In the latter case, the invariant chaotic set is the synchronization manifold. 
As the coupling strength varies over a given interval, the coupled map lattice presents 
UDV. The blowout bifurcation corresponds to the situation whereby the synchronization 
manifold loses transversal stability. A non-synchronized trajectory very close to this manifold 
experiences the same kind of kicks towards and away the synchronization manifold. The 3/2- 
exponent comes from the fact that the on-off intermittency comes from a system parameter 
being perturbed in the vicinity of a bifurcation point (in our case the blowout bifurcation). 
Our results provide further numerical evidence that the behavior of chaotic transients near 
a blowout bifurcation do present features universal to on-off bifurcation. 
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FIG. 7: Average inter-burst interval for models B and C with a = a* + <S, 6 = 10 -3 . Solid lines 
represent least-squares fit with slopes —1.492 and —1.480, respectively. 

The UDV-induced intermittency can also be viewed from the point of view of the shad- 
owability of trajectories in the chaotic invariant set Q. In practice, due to unavoidable 
factors as imperfect parameter identification and one-step truncation errors involved in the 
numerical procedures, a computer-generated trajectory thought to belong to Q will actually 
start off but very close to it. However, we know in advance that there exists a "true" chaotic 
trajectory in Q, which may or may not shadows the pseudo-trajectory obtained through 
numerics. 

The shadowing distance between the "true" chaotic trajectory at and the pseudo- 
trajectory initialized nearby is, at each instant, the pointwise distance between them in the 
phase plane. We usually compute the log-distance d n = log 10 z n (cf. Fig. 8) for the time- 
series of this variable for the same case depicted in Fig. 6. The existence of laminar intervals, 
for which the pseudo-trajectory is close to H, is equivalent to having a pseudo-trajectory 
which continuously shadows the "true" chaotic trajectory belonging to Q. Bursting is an 
observable manifestation of the lack of shadowability, while the lengths of the laminar in- 
tervals yield estimates for shadowing times. Hence, the properties of chaotic bursting are 
related to the statistics of shadowing distances and times. 
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FIG. 8: Evolution of the log-shadowing distances for the same time series depicted in Figs. 6(a) 
and (b). 

Figures 9(a) and (b) present numerically obtained statistical distributions of the shad- 
owing log-distances d n for models B and C, respectively, and different values of a, where we 
applied external kicks of level 5± = 10 ~ 16 as before. When a is far from the blowout bifur- 
cation (solid and dashed lines in Fig. 9) the (normalized) distribution height falls rapidly 
down to zero for shadowing distances less than 10~ 16 , as expected due to the noisy applied 
kicks, and decreases exponentially for higher shadowing distances V(d) ~ exp(— /cd), where 
k > 0. As a increases towards the blowout bifurcation, the slowdown rate k decreases until 
it becomes practically zero near the bifurcation (dotted lines in Fig. 9). This can be inter- 
preted from the fact that as the UDV effect is more intense near the blowout bifurcation, we 
have a progressive dominance of higher shadowing distances. This is in accordance with the 
greater content of transversely unstable periodic orbits as a is increased further from a*. 

The shadowing log-distances experience spikes of various heights, but remain in the im- 
mediate vicinity of the invariant subspace fi } until they burst chaotically and return to Vt. 
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FIG. 9: Probability distribution of log-shadowing distances for the map (a) B and (b) C, for 
different values of a and with addition of random kicks of strength 5± = 10 -16 along the transversal 
direction. 
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We define the shadowability time as the interval it takes for the pointwise shadowing dis- 
tance to grow until \z n \ > za = 0.001, which is a rather arbitrary choice but with little 
overall effect on our numerical results. The distribution of shadowing times is thus the same 
as the power-law distribution of laminar intervals P(r) depicted in Fig. 7. This fact can also 
be derived by integrating the probability distribution for shadowing log-distances V(d), in 
order to obtain the probability for a shadowing distance to be greater than i/a [7]. 

These probability distributions for the shadowing distances and times can be theoreti- 
cally justified from the statistical properties of finite-time Lyapunov exponents. A pseudo- 
trajectory starting off but near the invariant subspace will move either toward or away 
from Q for finite time intervals according to the corresponding time-ra transversal Lya- 
punov exponent: z^ n ~ z^ exip(n\_\_(n)). It follows that the log-shadowing distances satisfy 
dk^n ~ dk + n\±(n). When Q exhibits UDV, and especially when it is close to the blowout 
bifurcation, the time-ra exponents X±(n) act as random variables, what makes for an additive 
random process, with a diffusion rate given by the dispersion of the time-n exponents. Since, 
in general, the distribution of X±(n) is such that there is a different amount of positive and 
negative values, we obtain a biased random walk, with a drift towards a reflecting barrier, 
and which turns to be equal to the average transversal exponent. The external kicks we 
add are the reflecting barrier for the biased random walk. The distributions for shadowing 
distances and times may be thus obtained as equilibrium solutions of a diffusion equation 
taking into account these effects [7, 21]. 

VI. CONCLUSIONS 

We reported in this paper the parametric evolution of unstable dimension variability 
(UDV) in strongly non-hyperbolic dynamical systems of low dimensionality and having in- 
variant subspaces possessing a chaotic set. Two key situations are focused in this investiga- 
tion: the onset of UDV and its maximum intensity at the blowout bifurcation. In the former 
case, some unstable periodic orbit embedded in the chaotic set loses transversal stability, 
whereas in the latter, the chaotic set as a whole loses transversal stability. The continuous 
dependence of UDV on a bifurcation parameter can be described, in periodic-orbit theory, 
by the relative contribution of the unstable periodic orbits (UPOs) with a different num- 
ber of unstable directions. These contributions have been numerically tested through the 
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computation of the transverse finite-time Lyapunov exponents and a contrast measure that 
quantifies the relative contribution of the UPOs. 

The sign of the transverse finite-time exponents (which fluctuate about zero when UDV 
sets in) indicates that a trajectory of finite length is being pushed towards or away the 
invariant chaotic set on average. According to the periodic orbit theory, the abundance 
of positive exponents is related to the relative contribution of saddle and repeller UPOs 
embedded in the chaotic set. At the blowout bifurcation these contributions would match 
each other, verified for our models by using the computationally largest possible periodics 
for the UPOs. 

Moreover, in our analysis we deal with systems for which there may be UPOs with a 
larger number of unstable directions. For such systems, we found that the relation (22) 
between natural measure and the eigenvalues of UPOs is still valid, provided we stick to 
the original meaning, given in Ref. [3], for the expanding eigenvalues of the UPOs with 
unstable directions embedded in the invariant chaotic set, provided the latter is hyperbolic. 
In many systems, this is always the case (e.g. in Refs. [4] and [11]), and the difference in 
the adopted definition of expanding eigenvalue is immaterial. On the other hand, in the 
systems we deal with in this paper, exhibiting UDV for some parameter values, there may 
be subtle but important differences due to the possible existence of unstable directions off 
the chaotic invariant set. In any case, the agreement between periodic orbit theory and 
numerical experiments can be as good as one wishes, provided we choose UPOs with period 
large enough. 

We paid special attention to the mechanism underlying the onset of UDV which, in 
this work, differ qualitatively from that that holds for processes of the riddling bifurcation 
type [17] in two basic aspects. Firstly, for systems undergoing a transition via riddling 
bifurcation, the onset of UDV occurs through the loss of transversal stability of a low- 
period orbit embedded in the chaotic invariant set. In the mechanism we proposed in this 
paper, however, as we approach the onset of UDV the period of the orbits losing transversal 
stability goes to infinity. In the second place, in the riddling bifurcation scenario the onset 
of UDV involves the collision of three unstable periodic orbits: one transversely stable orbit 
embedded in the invariant chaotic set, and two other transversely unstable orbits outside 
the invariant set, such that there are unstable orbits existing before the onset of UDV. On 
the other hand, for the kind of transition shown in this paper, the onset of UDV is followed 
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by the creation of two transversely stable orbits after the onset of UDV. 

We have also considered the intermittent behavior of trajectories off but very close to 
the chaotic set. This intermittency is induced by UDV in the sense that the trajectory 
experiences bursting phases which steer them far away from the chaotic set whenever influ- 
enced by transversely unstable UPOs, but can return to the vicinity of the chaotic set under 
the chief influence of transversely stable UPOs. We obtained that, just after the blowout 
bifurcation, the probability distribution of inter-burst intervals scales according to the uni- 
versal 3/2 power-law characteristic of an on-off intermittency in the vicinity of the blowout 
bifurcation. These findings are according to previous results for UDV in coupled chaotic 
map lattices presenting transitions to non-synchronized behavior mediated by a blowout 
bifurcation [18, 31]. In terms of the shadowability of chaotic trajectories, we obtained an 
exponential distribution for both the log-shadowing distances and a power-law scaling for 
the corresponding shadowing times, in agreement with other systems studied and also with 
a theoretical model based on a biased random walk, with a reflecting barrier [7, 21]. 

To summarize, there are two key points to be emphasized in this paper. Firstly we 
reinterpreted the two-dimensional map studied in Ref. [11] in terms of UDV. While previous 
analyses focus on the blowout bifurcation, we also considered the onset of UDV, which occurs 
through a new mechanism: in the only case described in the literature the onset of UDV 
occurs through the loss of transversal stability of a fixed point [17], whereas in our case, as 
we approach the onset of UDV, the period of the orbits losing transversal stability tends to 
infinity (although this fact seems to depend on the chaotic dynamics in the invariant set). 
In the second place, we provided another higher-dimensional example which still allows for 
an analytical treatment for the counting of UPOs embedded in the invariant chaotic set. 
After describing for both systems the mechanism responsible for the onset of UDV and 
the point where it is the most intense, we considered the statistical behaviour of FTLE 
in order to show the observable consequences of the parametric evolution of UDV. The 
statistical quantities there obtained, the mean and variances, are useful to determine the 
time during which we expect shadowing trajectories to exist. Moreover, we found that the 
fraction of positive FTLE coincides, within the numerical accuracy, with the fraction of 
transversely unstable UPOs, what confirms the usefulness of the periodic orbit theory to 
describe UDV in a quantitative way. Although we have used low-dimensional model systems 
to draw these conclusions, essentially the same features are present in higher dimensional 
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dynamical systems. In particular, the systems we used presents the essential features of the 
synchronization (and transversal) dynamics of spatially extended systems like coupled map 
lattices or oscillator chains. 
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APPENDIX: Computing the first orbits to lose transversal stability 

We argued, in Section IV, that the onset of UDV occurs when the first UPO embedded 
in the chaotic invariant set loses transversal stability. In this Appendix we show, using a 
linear stability analysis, that the first orbit to lose transversal stability in model B has a 
diverging period and does so just after a = a c = 1.0. A similar analysis holds for model C 
and will not be presented here. 

Defining the product for all points of a period-p orbit x p (j) 

p-i 

Xpti) = II ( 2 %0') mod 1), (31) 
we can write the eigenvalue of this orbit along the transversal direction as 

e ± (J) = arx p {j), (32) 

such that this UPO will be transversely unstable (stable) if > 1 (< 1). Since < x < 1, 
it turns out that X p (j) is always less than unity and there are no transversely unstable UPOs 
when a < a c = 1. We are interested to find the period of the first orbit to lose transversal 
stability at a c = 1, hence we need to compute the value of C£fi rs t{p) f° r which the first 
period-p loses transversal stability and consider the value of p for which otf irst approaches 
a c = 1. 

We defined in Section IV, C£fi rs t(p) as the smallest value of a satisfying a p X p (j) > 1 for 
some period-p UPO, or 

ctfirM > (X' P )~ 1/P , (33) 
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where X f p = sup{X p (j)}. The orbit which initiates at the point x p (j = 2 P ~ 1 — 1 = J) is the 
only one possessing the remaining p — 1 points greater than 1/2, hence X f p = X P (J). There 
follows that (where each factor is to be understood as modi) 
'2P- 1 - 1\ (2 P - 2\ f2 pJrl - 2 2 \ f2 pJr2 - 2 3 



K 



2p - 1 J \2p - 1 

p / ™-i \ p 



n i 



•i=l 

in such a way that 



2 P ~ 



2 p - 1 



n i 



a 



first 



i=l 



(P)> 



2P-1 J V 2p - 1 
1 2 p 



2 2p- 



2 p- 



2^ - 1 



2 i 2p - 1 



n 



.i=l 



1 2 P 
2* 2? - 1 



(34) 
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Let di be the zth term of the product above. 



1 



1 2 p 



2 l 2 p - 1 

Since — > 1 when z (consequently p) goes to infinity, Eq. (36) furnishes a convergent result. 
Since the are always strictly positive and < a i+1 < 1, this product converges for some 
Kp such that < k p < 1 and Eq. (36) is rephrased as 



®first(p) > Kp 



(0 < K p < 1). 



(37) 



For large periods, since a p « 1, there follows that k- p+ i « and 

1) < OLf irs t{p), 

meaning that orbits with longer periods lose transversal stability earlier than orbits with 
shorter periods, as a increases. We remark that in other models, like that considered in Ref. 
[17], the orbits with lowest periods lose transversal stability first. In the p — > oo limit we 
find that 

lim a first (p) > 1, (38) 

p — *OG 

and the onset of UDV, at a = 1 + 5 with 5 —> 0, corresponds to the loss of transversal 
stability of an infinite period orbit. 
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